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ABSTRACT 

Microlensing event MACHO 97-SMC-l is one of the rare microlensing events 
for which the source is a variable star, simply because most variable stars are 
systematically eliminated from microlensing studies. Using observational data for 
this event, we show that the intrinsic variability of a microlensed star is a powerful 
tool to constrain the nature of the lens by breaking the degeneracy between the 
microlens parallax and the blended light. We also present a statistical test for 
discriminating the location of the lens based on the \ 2 contours of the vector 
A, the inverse of the projected velocity. We find that while SMC self lensing is 
somewhat favored over halo lensing, neither location can be ruled out with good 
confidence. 

Subject headings: gravitational lensing — stars: variables: other 



1. Introduction 

Microlensing parallax measurements are a potentially powerful way to constrain the 
nature of the lenses. For typical events, the only measured parameter that is related to the 
underlying physical properties of the lens is the Einstein timescale £e- However, this relation 
is rather indirect, 

a 

*E = — , OE = V^^rel (1) 

fJ'iel 

where 7r re i and // re i are the lens-source relative parallax and proper motion, $e is the angular 
Einstein radius, M is the mass of the lens, and k = 4G/c 2 AU ~ 8.14 mas Mq 1 . That is, £e 
is effectively a combination of three parameters describing the lens and the source: M, 7r re i, 
and fj, re \. If one can also measure the Einstein radius projected onto the observer plane, f E , 
or equivalently the microlens parallax, 
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then this three-fold degeneracy can be reduced by one dimension, allowing better constraints 
on the mass. Moreover, one can disentangle the mass from the kinematic variables and so 
extract a purely kinematic quantity, the projected velocity, 

f E AU/! r cl , oN 

V = T = ~Z • ( 3 ) 

This quantity is particularly useful for understanding the nature of the lenses detected toward 
the Magellanic Clouds (MCs), which is currently under debate. The MACHO collaboration 
(Alcock et al. 2000) detected 13 to 17 such events and argued that the majority of these 
were due to MACHOs making up about 20% of the Milky Way dark halo, whereas the EROS 
collaboration (Afonso et al. 2003; Tisserand & Milsztajn 2005) argued that their relative lack 
of such detections was consistent with all the events being due to stars in the Milky Way disk 
or the MCs themselves. All events detected so far toward the SMC have been very bright, 
even though the detection of fainter events was expected. We recognize this statistical issue 
but we will not address it any further on this work. One way to determine the nature of 
individual lenses is to measure their parallax 7Te and so their projected velocity v. Boutreux 
& Gould (1996) showed that this quantity differs greatly depending on the lens population: 
~ 50kms _1 for lenses in the Galactic disk, ~ 300 km s^ 1 for halo lenses, and ~ 2000kms~ 1 
for lenses within the MCs themselves. 

However, to date it has been possible to measure microlens parallaxes for only a handful 
of events (Poindexter et al. 2005 and references therein). Measurement requires that one 
compare fg to some "standard ruler" in the observer plane, which must therefore be of com- 
parable size to this quantity, i.e., of order 1 AU. While 7Te could be routinely measured by 
comparing the microlensing event as observed from the Earth and from a satellite in solar 
orbit (Refsdal 1966), the only standard ruler generally available for ground-based observa- 
tions is the Earth's orbit itself. Hence, unless the event takes a substantial fraction of a year 
(in practice t^ ^ 90 days), it is rarely possible to measure its parallax. 

The problem is that the microlens parallax is actually a 2-dimensional vector quantity, 
7Te, whose magnitude is tte and whose direction is that of the lens motion relative to the 
source. The component of 7Te parallel to the direction of the Sun at the peak of the event, 
7Te,||, induces a distortion in the light curve that is asymmetric with respect to the peak, 
while the other component, 7Te,±, induces a symmetric distortion. While 7Te ,ii is relatively 
easy to measure, even for comparatively short events (Gould et al. 1994), 7Te,j_ is subject to 
various degeneracies, both discrete and continuous (Smith et al. 2003; An et al. 2004; Gould 
2004). 

The fundamental problem is that changes in any one of four of the five parameters that 
describe standard microlensing also induce symmetric distortions in the light curve, and the 
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combination of these can often mimic the effects of 7Te,j_. The light curve of a standard 
(non-parallax) microlensing event is modeled by 

fk(t)=fs,kA[u(t)]+f h>k , (4) 

where the magnification A is given by (Einstein 1936; Paczyhski 1986), 

At ^ U 2 + 2 

u\]u l + 2 
and u is the source-lens separation in units of 8-e, which is given by the Pythagorean theorem, 

u = ^/{3 2 + r 2 (6) 

in terms of r and the impact parameter /3, 

T(t) = ^, P = tie. (7) 

Here fk(t) is the observed flux, / s ^ is the flux of the microlensed source, and /b,fc is the flux 
from any unlensed background light, all as observed at the fcth observatory. Inspection of 
these equations shows that changes in uq, £e, /s and /b all induce effects on the light curve 
that are symmetric with respect to to- 

If it were somehow possible to measure /b independent of the light curve, then the 
problems posed by these degeneracies would be greatly reduced. First, since the baseline 
flux, /base = /s + fb, is generally determined very precisely, measurement of fb immediately 
yields / s . Next, since the peak flux / pca k is also usually well measured, one can then also 
directly determine u from A(u ) = (/ pea k — fb)/ fs- This implies that only t E and 7r E ,j_ 
must really be simultaneously determined. Since changing these induces rather different 
light-curve distortions, it is relatively easy to disentangle them. 

To date, three methods have been proposed to determine fb- astrometry (Alard et al. 
1995; Ghosh et al. 2004), precise imaging (Han 1997), and manipulation of a series of images 
using image subtraction (Gould & An 2002). These methods can all help, but they are 
fundamentally unable to resolve blended light due to a binary companion of the source that 
does not participate in the microlensing event. Such companions could lie at separations of 
order 10 AU, which corresponds to angular separations of order 1 mas or smaller. For the 
faint sources typical of microlensing events, this is too close to be resolved using current or 
foreseen instruments. 

Here we propose a new method to resolve blended light. The method requires that the 
source be a regular variable, which is fairly rare for sources in typical microlensing fields. 
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However, the method does not suffer from the limitations of the other techniques: its ability 
to disentangle the source from the blended light works equally well regardless of their angular 
separation. We apply this technique to MACHO-97-SMC-1, which was a relatively long event 
whose source exhibited regular variability with an amplitude of about 3%. We show that this 
variability significantly constrains the blending and so improves the precision of the parallax 
measurement. 



2. Data 

MACHO 97-SMC-l was discovered by the MACHO collaboration (Alcock et al. 1997) 
and independently by the EROS collaboration (Palanque-Delabrouille et al. 1998). The 
OGLE collaboration also observed the event during its late decline in its OGLE-II phase 
(Udalski at al. 1997a) and at baseline in its OGLE-III phase (Udalski 2003). 

MACHO observations were mostly made at the dedicated 50" Great Melbourne telescope 
at Mount Stromlo, Australia from June 1993 to January 2000 with simultaneous imaging in 
"Vm" (4500-5900 A) and li R M " (5900-7800 A) passbands, but also at CTIO in the Johnson- 
Cousins R filter from May through November 1997. 

EROS observations were made at the dedicated lm MARLY telescope at the European 
Southern Observatory at La Silla, Chile from July 1996 to February 2003 with simultaneous 
imaging in U V E " (4200-7200 A, peak at 5600 A) and "7 E " (6200-9200 A, peak at 7600 A) 
wide passbands. 

OGLE observations were made at the 1.3 meter Warsaw telescope at Las Campanas 
Observatory, Chile in Cousins I band from June 1997 to November 2000 for the OGLE-II 
phase and from June 2001 to August 2005 for the OGLE-III phase. 

Because the OGLE-III observations began after the event was over, they cannot be used 
to directly constrain the microlensing event. We incorporate them only to better understand 
the variability. To this end, we fix the blending at a value similar to that of OGLE-II, since 
otherwise this parameter would be completely degenerate. 

Afonso et al. (1999) showed that the source star is an intrinsic variable with a period of 
5.126 days. They fit the event simultaneously for microlensing and intrinsic variations and 
used OGLE-II data to demonstrate that some of the apparent source in the EROS images 
was unrelated blended light and to better define the period of variability. The amplitude of 
these oscillations is so small (around 3% of the total flux) that the event was detected by 
EROS because it managed to survive their variable-star cuts. The MACHO collaboration 
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did not initially detect this event, mainly because their "main sequence variable" cut was 
still based on the LMC color-magnitude diagram and had not been corrected for the SMC. 
It was rather detected by coincidence while extracting variable stars for further study. 



3. Parametrization 

When a variable star is microlensed, the amplitude of the oscillations will be magnified 
by the same factor as the star's flux, completely independent of the amount of blended light. 
For example, if the star's variability period is much shorter than the event timescale, the 
ratio of the amplitude at baseline and at some magnified moment will yield a measurement of 
the magnification at that stage. Since this measurement is independent of the source's flux 
and of the blended light, knowing the magnification at some point allows one to solve the 
2x2 system of equations formed by equation (4) at these two moments, thus disentangling 
the source from the blended light and so eliminating the degeneracy between the parallax 
and the blending. 

In practice, measurements of the microlensing parameters are done by fitting a mi- 
crolensing model to the observations, eliminating the requirement for the variability period 
to be much shorter than the event timescale. In this process, the magnified variability am- 
plitude of the source will yield an extra constraint on the amount of blending to fit for the 
magnification at each stage, thereby breaking the forementioned degeneracy. 

In this paper we fit the event to a standard microlensing model (eqs. [4]- [7]) plus 
parallax and combined with a variability model for the intrinsic variations of the source star. 
Following the procedure outlined by An et al. (2002), we model the effects of parallax by 
modifying equation (7), 

P(t) = Uo + 6p, T(t) = t -^ + 6r, (8) 

where 

[5r(t),5f3(t)] = [7z E -As(t),7r E xAs(t)]. (9) 

and the parallax is evaluated in the geocentric frame. Here, As is the offset position of the 
Sun projected on the plane of the sky in the geocentric frame. The two components of 7Te 
are free parameters of the model that we fit to the data. As in equation (7), to is the time 
of the maximum, uq is the closest angular distance between the source and the lens in units 
of the angular Einstein radius, 8 E , and t E is the Einstein timescale of the event. 

Parallax models are generically subject to a four-fold discrete degeneracy (An et al. 
2004; Gould 2004; Poindexter et al. 2005). Two solutions are related by the "constant 



acceleration" degeneracy (Smith et al. 2003), which in the geocentric frame sends uq — ► 
—■u and changes other parameters relatively little. The two other solutions are related by 
the "jerk-parallax" degeneracy, in which the jerk of the Earth's motion can masquerade as 
a parallax effect (Gould 2004). Palanque-Delabrouille et al. (1998) found both constant- 
acceleration solutions, although they did not express them in exactly the terms presented 
here. 

To describe the intrinsic variation of the star, we will begin by following Palanque- 
Delabrouille et al. (1998), who modeled the variability as a sine function, leaving as free 
parameters its angular frequency Q, its phase 0, and its amplitude e k as measured at each 
observatory. Hence, the light curve is given by 

/*(*) = fs,k A(t) [1 + e k G(nt + </>)} + /„,* , G(x) = sinx. (10) 

Our approach is very similar, but with one minor modification. After we initially fit the light 
curve to equation (10), we fit the residuals to a fourth-order Fourier expansion to derive a 
new G(x), which differs modestly from a sine function (see Fig. 1). 

Udalski et al. (1997b) suggested that the light curve oscillations may be due to ellipsoidal 
variability of the source. If so, this would imply that the light curve should be subject to 
"xallarap" (binary-source motion) as well as parallax (Earth-motion) distortions. We will 
examine this possibility more closely in Appendix A, but for the present we simply treat 
these oscillations as an empirical fact of unknown origin. 



4. Results 

Using the model described in § 3, we fit simultaneously the data from all the observa- 
tory/filter combinations (7 in total). The x 2 of the best fit, in total and by observatory/filter, 
is shown in Table 1 for the positive and negative uq solutions. It should be noted that the 
errors are scaled (with scaling factors 1.92, 2.65, 1, 0.79, 1.1, 0.75 and 1 respectively for each 
observatory as listed in Table 1) so that the best fit x 2 is approximately equal to the number 
of data points in each observatory /filter. Also, the best fit parameters are listed in Table 2 
with their respective errors for both solutions. We searched for the jerk-parallax solutions 
according to the prescription of Park et al. (2004). We indeed found two additional minima 
on the x 2 surface, but these are excluded at the 12-0 and 28-0" levels respectively, and so 
they will not be further considered. 

Parallax parameters are of special interest to constrain the lens projected velocity and 
mass. Their best fit values show that there is very little, if any, parallax effect. In Figure 
2, we show x 2 contours for both (=bu ) solutions. Although the x 2 °f the best fits do not 
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differ significantly, the errors are a factor 2-3 smaller for positive Uq. In each case, the best fit 
variability period for the star is 5.1252±0.0001 days, compatible with, but much more precise 
than, the ones found by Palanque-Delabrouille et al. (1998) and Udalski et al. (1997b). 

Palanque-Delabrouille et al. (1998) raised the question of whether it is the microlensed 
source that is varying or it is the blended light. Their data alone permitted them to dis- 
criminate against the second option just at the 2.5-cr level. By including the data from all 7 
observatories we find that the variable blended-light model is ruled out at the 1-a level, so 
it will not be further considered. 

In order to evaluate the role of the intrinsic source variations in constraining the blended 
flux, and so 7Te,±, we remove the best-fit variation from the data and then refit these adjusted 
data using a standard microlensing model with parallax but with no oscillations. The x 2 an d 
the best fit parameters are summarized in Tables 1 and 2 respectively, for the ±Uq solutions. 
Note that the best fit for the parameters that are correlated with 7rE,_i_, namely uo, ts, f s ,k 
and /b,fc, have changed somewhat compared with the variability- model best fit, but the main 
effect is that their errors have increased by a factor of 2-3 while the x 2 °f the best fit is 
virtually the same. This is reflected in the difference between the left and right contour plots 
of Figure 2, which respectively do and do not take variability into account. The contours on 
the right are strongly stretched in the direction of tie,± (note that 7Teji points just 2°.2 north 
of west, so 71e,± is virtually the same as 7Te,n)- This demonstrates that if the source star is 
a regular variable, the degeneracy between 7Te,_l and /b can be broken and their values can 
be better constrained. 

Wyrzykowski et al. (2006) searched the entire OGLE-III database up through the 2004 
season for microlensing events with variable baselines. Out of about 1400 events, there were 
21 with periodic variable baselines and 111 that showed a baseline with irregular variability. 
The latter are not well suited for this method, in part because apparent "microlensing events" 
on irregular variables can in principle be just a manifestation of their variability and in part 
because the amplitude of variability during the event cannot be precisely determined from 
the baseline variability. Naively, one would expect that half of the periodic ones were due 
to a variable source and half due to variable blend, which translates into an expectation of 
a little less than 1 regular variable source microlensing event for every 100 events detected. 
When considering also the probability of detecting parallax asymmetries on the light curve, 
the number of events to which the method presented in this paper is applicable is somewhat 
smaller. So, even though events like MACHO 97-SMC-l are rare, there should be a fair 
number of cases for which the analysis described in this paper would prove useful. 
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5. Lens Location 

A better constraint on the parallax parameters also gives a better constraint on the lens 
location, since the projected velocity of the lens depends directly on them (eqs. [2] and [3]). 
Figure 3 shows the x 2 contours of the vector A, 

A ^ = ^TT' (11) 

whose direction is the same as the projected velocity but with the inverse of its magnitude. 
We choose this vector instead of the projected velocity because it is well behaved at high 
velocities and makes the contours clearer. Indeed if the best fit for t E were independent of 
7Te, the A contours would differ from the 7Te contours just by a scale factor. Note, however, 
that these x 2 contours are made in the frame of the Sun, rather than in the geocentric frame 
in which we conducted the parallax analysis. That is, we convert 

Vhel = V geo + V ffi> _L, (12) 

where v©,± = (29.75,3.10) kms -1 is the offset between the two frames. The inner circle 
shows v = 2000 kms -1 while the outer one shows v = 300 kms -1 . 

There are three possible locations for a lens that magnifies a star in the SMC: it is located 
either in the Milky Way disk, the Milky Way halo (halo lensing) or in the SMC itself (SMC 
self lensing). However, a disk location is immediately ruled out by the high projected velocity. 
At first sight, since the minimum x 2 is located around a projected velocity of 1000 kms -1 
for both uq solutions, SMC self lensing appears more probable, because this would seem to 
be a too high velocity for a halo lens. To quantitatively discriminate between these two 
possibilities, we compare the mean maximum likelihood parameter, (exp[— Ax 2 /2]), for halo 
and SMC lenses based on the x 2 contours of their A vectors and their respective velocity 
distributions. We consider SMC self lensing and halo lensing in turn. 



5.1. SMC Self Lensing 

We begin by assuming that SMC lenses are Gaussian-distributed, namely, that they 
have a projected velocity probability distribution given by 

P(v) dv x dv y oc e~( v/a) /2 v dv North dv Eas t, (13) 

where a is the velocity dispersion of the lenses. Then, the mean likelihood parameter is 

<e- A * 2/2 > = fe- A ^^e-^ 2 / 2 vd^ orth dv Eas J [ e^ 2 / 2 v dv mrth dv East . (14) 
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The typical projected velocity of a SMC lens is not very clear but should be of the order 
of v > 1000 kms" 1 , so we consider two different estimates for the velocity dispersion: o = 
1000/ V2 kms -1 and a = 3000/a/2 kms -1 . As mentioned above, the x 2 calculation is more 
stable in A space, so we evaluate equation (14) in this space: 

(e- A * 2/2 ) = je' A ^ A ^ e~^ 2 /* A" 5 rfA North rfA East / f e^/°? I* A ~ 5 rfA North dA East , 

(15) 
where A% 2 (A) is taken directly from the data of the contour plots of Figure 3. 



5.2. Halo Lensing 

The calculation of the mean likelihood parameter for halo lensing is very similar to 
the one performed for SMC self lensing in the previous section, but we must allow for two 
considerations that were implicitly neglected before: the spatial distribution of the lenses 
and the effect of the Sun's and the source's motions on the projected velocity of the lens. 
We assume that halo lenses are distributed in an isothermal sphere, 

p(r) = 4T^' (16) 

where r is the lens Galactocentric distance and where we adopt a = 5 kpc for the core 
radius and assume that the Sun is at Rq = 7.6 kpc. The velocity probability distribution 
is then weighted by the density of lenses at each distance and by the size of its respective 
Einstein ring, and integrated over all possible distances, giving a mean likelihood parameter 
of 

( e - A x 2 / 2 ) = J exp [-A X 2 ( A)/2] g(A, D L ) d D L dA North dA East / f g{A, D L ) d D L dA NoTth dA EasU 

(17) 
where 

\2 



9(A , Dl) = °xpHW t »V2] A- D „ P (D L) 1°*" 



2Tr{a(D L )Y- M " "' V D s ' v ' 

The velocity dispersion is now given by 

HD L) = ^ ^. (19) 

with v rot = 220 kms -1 . Here, Di and Ds = 60 kpc are the distances to the lens and the 
source, while D^s = Ds — D^. The extra a' 2 factor in equation (18) is due to the fact that 
we cannot ignore the normalization of each Gaussian in this case. Note that in equation 
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(18), vl is regarded as an implicit function of A, which must be made explicit in order to 
quantitatively evaluate equation (17). Also note that there is an extra factor of D^s on 
equation (18), since the probability should be weighted by v and not by v, as was done in 
the previous section. 

The second consideration we must take into account for halo lensing is the effect of the 
motion of the Sun and the source on the projected velocity of the lens. In equation (18) we 
are assuming that the velocity distribution of the lenses is isotropic in the Galactic frame 
but we must also include the velocity offset induced by the motion of the Sun and the SMC 
with respect to the Galaxy. The three dimensional velocity of the Sun in UVW coordinates 
in the Galactic frame was measured by Hogg et al. (2005) to be 

v Q = (10.1, 224, 6.7) kms" 1 , (20) 

where U points from the Sun to the Galactic center, V points towards the Sun's Galactic 
rotation and W points towards the Galactic north pole. The velocity of the SMC is not 
well measured, but since it is most probably orbiting the LMC, it should have a very similar 
three dimensional velocity. Hence, we adopt the three dimensional velocity of the LMC for 
the SMC, which was measured by Kallivayalil et al. (2005) to be 

«smc = (-86, -268, 252) kms" 1 , (21) 

in the same coordinates system. Projected on celestial north and east coordinates, these 
velocities are 

(t> 0)N , w 0)E ) = (119, -133) (22) 

(vsmc,n, «smc,e) = (-296, 229). (23) 

The projected velocity of the lens relative to the observer-source line of sight is then 

v = v L - v (D L ), (24) 

with 

v (D L ) = v Q + —^v SMC . (25) 

Solving equation (24) for vl and substituting into equation (18), we find that for the positive 
Mo solution, the ratios of the mean likelihood parameters of the SMC to the halo are 7 and 3 
for the larger and smaller velocity dispersion respectively, while for the negative uq solution 
they are 3 and 2 respectively. This means that while SMC self lensing is somewhat favored, 
neither of the locations can be ruled out with good confidence. 

In the left panel of Figure 4 we show the A probability distribution for halo lenses 
magnifying stars in the SMC, integrated over all D^. There are two probability peaks with 
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similar A (or speed) but roughly opposite direction. This is due to the velocity offset, which 
flips sign as Dl increases. Also, for completeness of the statistical method proposed in this 
paper, we show in the right panel of Figure 4 the same probability distribution but for halo 
lenses causing microlensing events toward the LMC, in which an asymmetric probability 
peak can also be seen. 

It should be noted that throughout this statistical test, we have assumed that there is 
equal prior probability that a microlensing event is caused either by self lensing or by halo 
lensing. There is no evidence for assuming something different, so the ratio of the mean 
likelihood parameters is actually the ratio of the probabilities of having this event caused by 
self lensing and by halo lensing. 



6. EROS Fixed Blending 

OGLE observations showed the presence of an optical companion to the source of MA- 
CHO 97-SMC-l that was unresolved by EROS. Udalski et al. (1997b) estimated the light 
contribution of the companion to be 23%-28% of the total / band flux. This compares with 
28 ± 12 % for the fits shown in Table 2, which tends to confirm that the OGLE-resolved 
source is the true microlensed source (although it remains possible in principle that there is 
a small component of blended light within this source). In this section we therefore assume 
that the blended flux in the EROS Ie band is exactly 24% and repeat the analysis of § 4. 
Such a strong constraint should remove the degeneracy between /b and 7Te,±, perhaps leading 
to a better determination of the lens location. 

Refitting the model described in § 3 under this assumption, we obtain the x 2 contours 
of the parallax parameters for both (±«o) solutions shown on Figure 5. Both sets of contours 
are significantly less elongated in the 7Te,_l direction when compared to their analogs in Figure 
2, showing that the degeneracy is mostly broken, as expected. 

A better constraint on the microlensing parallax allows for a better constraint on the 
projected velocity. Following the analysis detailed in § 5, we show the x 2 contours for the 
components of vector A in Figure 6. These contours are also strongly compressed, so a better 
determination of the location of the lens might be possible. 

Employing the statistical test detailed in § 5, we find that, for the positive Uq solution, 
SMC self lensing is favored by factors of 27 and 33 for the larger and smaller SMC assumed 
velocity dispersions respectively. However, for the negative uq solution we find that the 
factors are 1 and 3, respectively. This result is very important, since it shows that even 
if we impose a strong constraint on the blended light, the location of the lens cannot be 



-14- 



determined with good confidence for this event. 



7. Discussion 

Variable stars are systematically eliminated from microlensing studies, but if one of 
them undergoes a microlensing event, its variability might actually be useful to constrain 
the event's nature. 

One such an event, MACHO 97-SMC-l, passed the variable-star cuts used by the MA- 
CHO and EROS collaborations. Using several data sets of this event taken by the MACHO, 
EROS and OGLE collaborations, we have shown that if the magnified source star is a regular 
variable, the parallax parameter in the direction perpendicular to the projected position of 
the Sun, 7Te,_l, can be constrained much better than if it were not a variable: the variability 
breaks the degeneracy between 7Te,± and the blended light. Figure 2 shows that, if we remove 
the variability from both the model and the data, the width of the 7Te contours are almost 
doubled in the 7Te,_l direction while they remain constant in the tte,\\ direction. 

We have proposed a statistical test to discriminate between the possible lens locations 
(halo or SMC) based on their velocity distributions and the \ 2 contours of A, the inverse of 
the projected velocity. For MACHO 97-SMC-l, the test revealed that the lens is more likely 
to be located in the SMC itself rather than in the Galactic halo but that it is not possible to 
rule out the latter location. This test might be useful to discriminate quantitatively between 
possible lens locations of other microlensing events for which a parallax measurement can 
be performed. 

We also showed that, even imposing a strong constraint on the flux, indicated by OGLE 
observations, the location of the lens still cannot be determined with good confidence. 
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97-SMC-l source star. Work by RJA and AG was supported by grant AST-0452758 from 
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A. MACHO 97-SMC-l as an Ellipsoidal Variable 

Udalski et al. (1997b) suggested that the source star is an ellipsoidal binary system, 
based on its nearly sinusoidal variations. In such systems, the binary separation is very 
small, and either one or both components are ellipsoidally distorted. Due to the interactions 
between them, their rotation and orbital periods are usually synchronized and the orbits are 
most likely circularized (see Beech 1985 and references therein). The star appears brighter at 
quadrature than at conjunction for two reasons. First, it has a larger surface area facing the 
observer. Second, the mean surface gravity of the exposed area is higher, implying that the 
surface is hotter and thus has a higher surface brightness. This leads to sinusoidal variations 
with a period equal to half of the orbital period. 

Sahu & Sahu (1998) showed that, if MACHO 97-SMC-l is a binary system, it is a single- 
lined spectroscopic one. As outlined by Morris (1985) and reviewed and updated by Morris 
& Naftilan (1993), it is possible to estimate the mass of the secondary star as well as its 
radius and other orbital parameters for a single-lined spectroscopic ellipsoidal binary using its 
light curve, the mass of the primary and its radial velocities (for double-lined spectroscopic 
binaries the mass of the primary can also be derived). The companion star should be dim 
enough to not be detected in the spectrum but massive enough to produce an ellipsoidal 
distortion on the primary. 

To relate the amplitude of ellipsoidal variations to the binary parameters, we use equa- 
tion (6) of Morris (1985), 

qRlsm 2 i 3.070 AM 1 (3 - «i) 



A 3 (n + l)(15 + ui) 



(Al) 



where q is the binary mass ratio (1712/1711), i is the orbital inclination, R\ is the mean radius of 
the primary, A is the semi-major axis of the orbit, AMi is the mean peak-to-peak amplitude 
of the primary star alone, U\ is the linear limb-darkening coefficient for the primary and T\ 
is its gravity darkening coefficient. These last two parameters (t\ and u\) depend on the 
wavelength of the observation, A. In our specific evaluation we will assume a Ve filter, but 
note that the results are essentially independent of this assumption. Equation (6) of Morris 
(1985) is accurate to first order and will suffice for this analysis. A more exact version of 
this equation is given by Morris & Naftilan (1993). 

From its position on the color-magnitude diagram (Fig. 7) the source is a late B main 
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sequence star (see also Palanque-Delabrouille et al. 1998; Sahu & Sahu 1998) with an ap- 
proximate mass of 5 M & . Using the empirical calibration of van Belle (1999) between the 
angular diameter and the (V — K) color yields Ri ~ 5R & if we assume that the value of 
(V — K) is twice the value of (V — I) = —0.05 and that -Dsmc — 60 kpc. Theoretical models 
predict a 15% smaller radius, but this will not be important for the results of this analysis. 
The orbital period P should be twice the variability period, so P = 10.2 days according to 
the variability period found in § 4. 

According to the tables of Al-Naimiy (1978), for the Ve filter u\ rs 0.35, while T\ is 
given by equation (10) of Morris (1985) 

= 1.43879 x 1Q8/AT! 

1 P l-exp(-1.43879x lOVATi) l ' 

where Ti is the effective temperature of the primary star (14000K in this case, Sahu & Sahu 
1998), A is in A and (5 = 0.25 for early type stars (von Zeipel's Law, von Zeipel 1924), and 
yields T\ = 0.55. Using Keplers Third Law to eliminate A in equation (Al) and solving for 
AMi, we find that 

AM, = 0.0094 sin 2 , f^-V ( ^-Y ( P V' ^ (A3) 

\5R Q J \5M Q J V 10 - 2da y s / m l +m 2 v ; 

From Table 2, AM ± w [5/ln(10)] e Vs = 0.063 mag for the U V E " filter. If we adopt the 
most extreme mass parameters, m-i ^> m\ and % = 90°, the right-hand side must increase by 
a factor of 7 from our estimate to reproduce the observed AM\. 

This type of binary system should also present xallarap distortions on its microlensing 
events light curve. In order to test this, we fitted the data with a standard microlensing 
model with xallarap and variability but no parallax and found out that, to the 3-<r level, 
Xe < 0.015, where xe is the xallarap parameter, 

Xe = -• (A4) 

That is, xe is the semi-major axis of the primary component, ai, in units of the Einstein 
ring radius projected in the plane of the source, fg. This small xallarap tends to reinforce 
the idea that MACHO 97-SMC-l is not a binary system and is therefore not an ellipsoidal 
variable. However, the small xallarap might simply be caused by a very large Einstein ring 
projected on the source plane, which cannot currently be ruled out since the location of the 
lens cannot be determined with good confidence. 
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Observatory 


Filter 


u > 


u < 


u > 

NoVar 


u < 
NoVar 


MACHO 


V M 


306.6 


306.5 


305.6 


305.9 




Rm 


1057.7 


1056.7 


1052.0 


1053.1 


EROS 


V e 


724.5 


724.5 


724.5 


724.6 




Ie 


841.0 


841.9 


843.8 


844.1 


MACHO CTIO 


R 


74.9 


75.0 


75.9 


75.6 


OGLE II 


I 


315.9 


315.9 


316.4 


316.2 


OGLE III 


I 


367.3 


367.3 


367.3 


367.3 


Total 




3688.0 


3687.8 


3685.5 


3686.9 



Table 1: x 2 f° r each filter (for each observatory) for the different models and solutions. NoVar 
indicates the case for which the best fitted variability of the star has been taken out from 
the data. Note that, because of round-off, in some cases the \ 2 of the filters/observatories 
do not add exactly to the total. 
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u > 





u < 





u > NoVar 


u < NoVar 




Best Fit 


Error 


Best Fit 


Error 


Best Fit 


Error 


Best Fit 


Error 


t (days) 


1460.041 


0.338 


1459.936 


0.385 


1460.014 


0.329 


1459.910 


0.428 


u 


0.419 


0.048 


-0.423 


0.046 


0.239 


0.071 


-0.344 


0.077 


t E (days) 


124.299 


9.780 


129.059 


13.689 


188.727 


44.494 


158.647 


37.382 


7TE,N 


-0.017 


0.011 


-0.028 


0.020 


-0.057 


0.015 


-0.059 


0.032 


7TE,E 


-0.002 


0.003 


-0.005 


0.009 


-0.004 


0.003 


-0.005 


0.008 


SI (rad/days) 


1.226 


0.000 


1.226 


0.000 










t ,e (days) 


2.539 


0.025 


2.540 


0.025 


















MACHO 








fs,V M 


0.550 


0.090 


0.557 


0.088 


0.259 


0.096 


0.417 


0.128 


fb,V M 


0.253 


0.090 


0.246 


0.088 


0.543 


0.096 


0.385 


0.127 


e v M 


0.029 


0.004 


0.028 


0.004 










Js,Rm 


0.861 


0.141 


0.871 


0.138 


0.405 


0.151 


0.653 


0.200 


Jb^R^i 


0.360 


0.141 


0.350 


0.137 


0.815 


0.150 


0.568 


0.200 


€Rm 


0.023 


0.003 


0.023 


0.003 


















EROS 








fs,V E 


0.532 


0.087 


0.539 


0.085 


0.251 


0.093 


0.403 


0.124 


fb,V E 


0.236 


0.087 


0.230 


0.085 


0.518 


0.093 


0.365 


0.124 


e v E 


0.029 


0.005 


0.029 


0.005 










JsJe 


1.226 


0.201 


1.241 


0.196 


0.577 


0.215 


0.929 


0.285 


fb,I E 


0.487 


0.201 


0.473 


0.196 


1.136 


0.215 


0.784 


0.285 


eis 


0.024 


0.004 


0.024 


0.004 


















MACHO CTIO 








fs,R 


24.133 


4.126 


24.426 


3.950 


11.384 


4.134 


18.351 


5.581 


fb,R 


7.277 


4.155 


6.965 


3.967 


20.026 


4.101 


13.032 


5.574 


£i? 


0.030 


0.005 


0.029 


0.005 


















OGLE II 








fs,I 


0.734 


0.128 


0.737 


0.127 


0.347 


0.134 


0.551 


0.174 


fb,I 


0.204 


0.129 


0.201 


0.127 


0.590 


0.133 


0.386 


0.174 


<=1 


0.031 


0.006 


0.031 


0.006 


















OGLE III 








fs,I 


0.684 


0.001 


0.684 


0.001 


0.684 


0.001 


0.684 


0.001 


fb,I 


0.254 




0.254 




0.254 




0.254 




e/ 


0.033 


0.003 


0.033 


0.003 











Table 2: Best fit parameters with their respective errors for the different models and solutions. 
NoVar indicates the case for which the best fitted variability of the star has been taken out 
from the data. t 0e is defined as <f>/£l, where <p is the phase defined in equation (10). The 
time of the maximum, to, is measured in JD — 2449000, while the fluxes are measured 
in arbitrary units. The seemingly excessive number of decimal places included for each 
parameter are potentially important to reproduce these results, because the correlations 
among some parameters are high. 
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Phase 



Fig. 1. — Variability light curve of MACHO 97-SMC-l, derived from the residuals of the 
fit of equation (10) with G(x) = sin(x). The solid curve shows the variability model used 
in this analysis, which correspond to a fourth-order Fourier expansion fit of the residuals, 
while the dashed line shows a sine function, the model used by Palanque-Delabrouille et al. 

(1998). 
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Fig. 2. — Ax 2 contours (1,4,9,16,25,36,49) of the parallax parameters for positive (top) and 
negative (bottom) uq. The plots on the left show these contours for the model described in 
§ 3, a standard microlensing model with parallax and taking into account the variability of 
the source, while the ones on the right shows the parallax x 2 contours for the same model 
but with the variability removed from both the data and the model. The contours on the 
right are strongly stretched in the direction of 7Te,_l (which coincides almost exactly with 
tte,n) compared to the ones on the left, showing that the variability of the source breaks the 
degeneracy between 7Te,± and /b,fc, allowing to constrain much better the parallax. 
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Fig. 3. — Ax 2 contour (1,4,9,16,25,36,49) plots for A, the inverse of the projected velocity 
(in kms -1 ), for both ±u solutions. The small circle in the middle shows a projected velocity 
(in the frame of the Sun) of 2000 kms -1 while the other circle shows 300 kms' 1 , typical 
projected velocities for SMC and halo lenses, respectively. 



:/.: 



O 
O 

o 



X, 
o 

< 




6 4 2 -2 -4 -6 

A F (1000 kin/s)- 1 




2 0-2-4 -6 

A East (1000 kin/s)- 1 



Fig. 4. — Inverse projected velocity A probability distribution of halo lenses causing mi- 
crolensing events toward the SMC (left) and the LMC (right). The velocity offset discussed 
in § 5.2 makes the distributions highly asymmetric, contrary to the case for SMC self lensing. 
The contour levels are offset from one another by a factor of 5, with order, highest to lowest, 
bold solid, bold dashed, thin solid and thin dashed (of which there are several). 
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Fig. 5. — Ax 2 contours (1,4,9,16,25,36,49) for the parallax parameters assuming that 24% of 
the light detected in the EROS Ie band is due to blending. Compared to the left contours in 
Figure 2, they are strongly compressed in the 71e,i. direction, showing that by this assumption 
the degeneracy between /b and 7Te,_l is almost completely eliminated. 
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Fig. 6.— Ax 2 contours (1,4,9,16,25,36,49) for A assuming that 24% of the light in the EROS 
Ie band is due to a blend. As discussed in § 6, even though the contours are now nearly 
circular, the location of the lens cannot be determined with good confidence. 
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Fig. 7.— Color Magnitude Diagram of the OGLE-II SMC field centered on MACHO 97- 
SMC-1. The star marks the position of the source of the MACHO 97-SMC-l event. This 
indicates that it is a late B main sequence star, compatible with the spectroscopic determi- 
nation of Sahu & Sahu (1998). 



